On the electron-energy loss spectra and plasmon resonance in cuprates 
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The consequences of the non-Drude charge response in the normal state of cuprates and the effect 
of the layered structure on electron-energy loss spectra are investigated, both for experiments in the 
transmission and the reflection mode. It is shown that in the intermediate doping regime the plasmon 
resonance has to be nearly critically damped as a result of the anomalous frequency dependence of 
the relaxation rate. This also implies an unusual low-energy dependence of the loss function. Both 
facts are consistent with experiments in cuprates. Our study based on the t-J model shows good 
agreement with measured plasmon frequencies. 



PACS numbers: 71.27.+a, 79.60.-i, 71.20.-b 

In metals the investigation of the plasmon resonance is 
used as a standard tool to extract the information about 
charge carriers. As several other electronic properties 
also the plasmon response is quite unusual in cuprates. 
This can be realized from the dielectric measurements 
jl| which reveal a very broad plasmon resonance in the 
electron-energy loss (EEL) function and hence a poorly 
defined plasma frequency ui p . One of the experimen- 
tal approaches to measure the plasmon dispersion w p (q) 
is the electron-energy loss spectroscopy (EELS) [||. So 
far the EELS has not been employed very frequently in 
cuprates, also due to the ambiguities in the interpreta- 
tion. The EELS in the transmission mode has been per- 
formed mostly for Bi 2 Sr 2 CaCu 2 8 (BSCCO) §{§, show- 
ing an optic plasmon with a quadratic dispersion. Sim- 
ilar spectra, restricted to q — 0, have been obtained for 
La2- 2; Sr 2 ;Cu04 (LSCO) by analysing the dielectric reflec- 
tivity H . Recently higher- resolution EELS spectra have 
been measured || for BSCCO in the reflection mode, 
displaying an acoustic-like plasmon at small wavevectors 
q. In both experiments the plasmon resonance appears 
heavily damped, although not overdamped. 

So far there have been only few theoretical studies of 
EEL, i.e. relevant to the quite specific electronic struc- 
ture of cuprates. It is has been realized early, in analogy 
with the layered electron gas Q] , that weakly coupled lay- 
ers of mobile electrons in CuO planar structures can lead 
in addition to optic plasmons also to acoustic plasmon 
branches || . Recently, an important relation of the EEL 
function and the related dielectric function e(q, uj) to the 
mechanism of the superconductivity in cuprates has also 
been stressed ||. 

Another important aspect follows from the observa- 
tion that the charge dynamics in the normal state of 
cuprates is anomalous and not consistent with the Drude- 
type relaxation. Charge as well as spin fluctuations have 
been instead phenomenologically described within the 
marginal- Fermi-liquid (MFL) scenario pC|Jllj |. In partic- 
ular, the latter has been shown to apply in the optimum- 



doping regime to the optical conductivity <j{uj), which 
can be expressed within a generalized Drude form with an 
effective (MFL-type) relaxation rate r _1 ~ 2ttX(uj + £T) 
pi} . More recently, the same anomalous dynamical be- 
havior has been found in numerical studies within the t-J 
model at intermediate doping 12-14)1, revealing a more 
precise form of this novel diffusion behavior. Closely 
related to the dielectric function e(q, uj) and the EELS 
in cuprates are calculations of the density fluctuations 
N(q,ui), studied within the t-J model both analytically 
|l5| and numerically 

The aim of this contribution is to discuss the conse- 
quences of the established non-Drude planar optical con- 
ductivity for the loss function and EELS experiments, 
both in the transmission and in the reflection mode. In 
particular we shall investigate the intrinsic damping of 
the plasmon mode. 

The microscopic models, such as the Hubbard model 
and the t-J model | p7[ , studied in connection with 
electronic properties of cuprates and other materials 
with strongly correlated electrons do not incorporate 
the long-range Coulomb interaction. The latter is nec- 
essary for the appearance of plasmon oscillations. It 
is quite straightforward to include the long-range inter- 
action within the random-phase approximation (RPA), 
where the dielectric function e(q, uj) and the dynamical 
density susceptibility x(q, uj) can be represented as 



e(q,w) 

x(q,w) 



i + V" q xo(q,w), 
i + Vqxo(q,^)' 



(i) 



where Xo(q7 w ) is the electron density suceptibility which 
includes the short-range correlations. In cuprates the 
hopping-matrix element for electrons between layers is 
rather small, as e.g. manifested in the large resistivity 
anisotropy, so that xo(q>w) should depend only on the 
planar component q||, i.e. xo(qi|>w). On contrary, the 
Coulomb matrix element V'q still depends on the 3D wave 
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vector q = (c\\\,q z ). The general expression for Vq has 
been given in Refs. |?],||]. For simplicity we restrict our 
discussion to small q, i.e. q z d < 1, qua < 1, where d, a 
denote the interlayer distance and the planar cell dimen- 
sion, respectively. In this regime one gets 



Vq = e 2 /e (e l \qf l +e±q%), 



(2) 



where en and e± are planar and interlayer u — > oo di- 
electric constants, respectively, due to screening of non- 
conduction electrons. The loss function is given by 0j 



Tm 



Vqlm X (q,oj). 



(3) 



In the following we shall focus on the long-wavelength 
limit q — > and on the regime of higher frequencies u ~ 



Here it follows from the continuity equation, 



Xo(q||,w) 



iq\a(u) 



(4) 



where <j(u) is the (complex) planar conductivity, which 
includes short-range correlation effects. 

The optical conductivity cr{u) — Recf(w) in cuprates 
has been studied extensively jlj , revealing anomalous be- 
havior consistent with the MFL scenario [ p~0|JXT| ] . There 
have been also numerous theoretical calculations of cr(u) 
at T = 0, using exact diagonalization of finite-size models 
for strongly correlated electrons in cuprates (ll| . Here we 
restrict our attention to the t-J model jl7|], which incor- 
porates the interplay between antiferromagnetic (AFM) 
spin correlations governed by the exchange J and the 
motion of electrons governed by the hopping parame- 
ter t, and the exclusion of double occupation of sites. 
Note that within this model the only relevant parameters 
are J/ 1 and the concentration Ch of holes introduced by 
doping. It is expected that the t-J model can represent 
(even quantitatively) quite realistic low-energy properties 
of electrons in cuprates [Q . Similar results should also 
follow for the Hubbard model, provided that the strong 
correlation regime J -C t, U <C t is considered. For a 
number of observables it is established that the inclu- 
sion of the next-nearest neighbor hopping t' can improve 
the agreement with experiments p"8] |. However there are 
indications that the qualitative behavior of cr(u) is less 
sensitive to the introduction of t' |l4j , therefore we omit 
such terms in this study. 

Recent numerical studies using the finite-temperature 
Lanczos method (for a review sec |Q) of the t-J model 
with realistic parameters J/t = 0.3 reveal results for 
cr(u), which are quite consistent with experiments in the 
normal state of cuprates. In particular, it has been found 
in the intermediate doping regime, 0.1 < < 0.3, that 
in a broad frequency range u < lo* ~ 3t a (to) follows a 
novel universal law [I13l , 



a{u) = C 



1 



UJ < UJ . 



(5) 



This anomalous diffusion has its origin in the fast decay 
of current-current correlations within a disordered spin 
background and seems to be generic for correlated sys- 



tems 1 19 , at least at intermediate temperatures T > T*. 

On the other hand, in a tight-binding model one can 
always express the complex conductivity ct(cj) in terms 
of a memory function M(u>) as 



ie 2 Kt/d 
J + M(uj) ' 



K 



2Nt' 



(6) 



where /C is the dimensionless kinetic energy density (per 
cell) at given T. For the case of a weak scattering 
|M(w p )| Up, Eqs.(0,§,^g) yield the plasma frequency 



eoeiicf 



(7) 



Ekin — (Hkin) is quite well established at T = from 
small-system studies Q , and it is expected to vary only 
slightly at low T < J @. Note that at very low 
doping, e.g. c/, < 0.05, assuming independent holes 
one expects Ekm oc c^. For higher doping the be- 
havior remains qualitatively similar, although the slope 
d\Ekin\/ dch is reduced. Taking into account our numer- 
ical values for Ekin obtained for J = 0.3 t and systems 
with N = 16, 18, 20 sites as well as the corresponding 
number of holes < Nh < 5, we get the results for 
u p (ch = Nh/N) shown in Fig. 1. In the evaluation of 
Eq.(0) we assume t ~ 0.4 eV and use parameters 
relevant for LSCO, i.e. en = 4 and d = 0.66 nm 



> 




FIG. 1. Plasma frequency uo p calculated within the t-J 
model assuming parameters relevant for LSCO: t = 0.4 eV, 
J = 0.3 t, too = 4 and d = 0.66 nm. 

It is evident that Eq.([|) cannot be modelled with small 
M(u). It is useful to present a(u) in the generalized 
Drudc form [til, 



a(u) 



r)(u)[u + ij(u)]' 



(8) 
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with an effective relaxation rate (inverse relaxation time) 
7(0;) and the mass renormalization f](u>), 



7(w) 



M"{lo) 

Tl(u) 



r){cj) = 1 + 



Af'(w) 



(9) 



Since r/(u> p ) ^ 1, it would be more appropriate to de- 
fine to* = Up/ J r](u>*), which should more directly ap- 
ply to the q z = plasmon resonance, relevant for op- 
tical mesurements and for the transmission EELS. Our 
numerical T > results for <j{lo) and consequently for 
M(u>) jl3],[l4} reveal however quite consistently for dop- 
ing concentrations 0.05 < ch < 0.3 that rj(u) p ) < 1, hence 

Up ~ «V 

A peculiar feature of the charge dynamics in cuprates 
is the non-Drude behavior of ct(lo), as e.g. represented 
by the MFL scenario , or more specifically by Eq. (||) , 
which should be valid at least near optimal doping. 
Within the broader MFL concept the dependence j(lo,T) 
0,0 is given by, 



7 = A(co + £T), 



A = 2ttA. 



(10) 



This behavior is obeyed in cuprates up to lo lo p . The 
dimensionless constant A obtained from numerical studies 
is not small. Estimates for various compounds fall in the 
range 0.5 < A < 0.9 @. 

Although A appears as a free parameter within the 
MFL proposal, this is not the case with Eq.(||), which 
also exhibits the MFL variation, Eq.(10), apart from log- 
arithimic corrections. Assuming that a {to > lo*) = 0, one 
can analytically evaluate <t(cj) and consequently M(lo). 
It is easy to express 7 in two regimes: 
a) For lo <C T <C lo* the expansion in lo yields 



A 



3 + 21nT/a 



b) while for T < w < w* we get 



2(l + lnw/T)' 



(11) 



(12) 



For the relevant range of T in the normal state of cuprates 
and the experimental range of lo (with reliable results), 
Eq.(|l2"|) should apply to experiments. The best overall 
fit of the numerical results (l2|,|l4j for T <S^t and u < t is 
achieved with A ~ 0.6 and £ ~ 2.7, also quite consistent 
with experiments 

If we insert the generalized Drude expression (g) and 
Eq.(|), into Eqs.ffl), we get 



e(q,w) = 1 



(13) 



where w p (q) is the effective plasmon frequency for a lay- 
ered system, 



a)p(q) 



(14) 



Note that for qu < q z the plasmon shows an acoustic 
dispersion, i.e. oj p — LL> p q\\e\\/q z ej_ oc q\\ [0J|]. 

With Eq.@ the EEL function (§ can be written as 



J( W ) = 



{rjLO 2 



- uj 2 ) 2 + j 2 rj 2 oj 2 



(15) 



Let us restrict our attention again to the intermediate 
doping regime where relations (p|,|l0|) apply. In the range 
to > T, mostly relevant to experiments and for the study 
of the plasma resonance in particular, we insert 7 ~ Xlo 
and assume rj(uj) ~ 1. Note that these simplifications 
could be possibly violated in the extreme acoustic limit 
Lo p (q) — * 0. The result is 



(X 2 - l) 2 + X 2 X 4 



LO 
LO n 



(16) 



It is evident from expression (16) that the width of the 
plasmon resonance is entirely determined by A, i.e. the 
resonance width is Alj/lo p ~ A. A further characteristic 
feature is the anomalous variation below the resonance, 
i.e. I (to < lj p ) oc lo 2 . 




FIG. 2. Loss function I (Eq. [16]) vs. x — uo/ui p for 
three different A-values, i.e. A=0.5, 1.0 and 2.0. The inset 
shows the ratio F/xo (solid curve) between the half- width 
r at half-maximum measured at the low-energy side of the 
peak and the renormalized plasma frequency xo = Lu max /tOp 
(dashed curve). 

Figure 2 presents I(u>) for three different values of A. 
For a typical value A ~ 1 this implies that the plasmon 
resonance is necessarily very broad, i.e. the halfwidth of 
the peak of I{lo) is of the order of the plasma frequency 
itself. This seems to be a new aspect of the anomalous 
behavior of strongly correlated electrons in doped AFM 
and cuprates. In addition, there is a noticeable downward 
shift of the peak position with respect to lo p . 
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Let us turn to the discussion of experimental results of 
cuprates. The loss function /(q, ui) has been determined 
directly by EELS measurements for BSCCO ||f| and for 
YBa2Cu3C>7 |2(J]. The characteristic asymmetric shape 
of I(u) and the strong damping of the plasmon peak, 
emerging from our analysis, appears to be consistent with 
these experiments 

The determination of the loss function I(to) from opti- 
cal q — > data turns out to be also difficult, in particular 
because of the limited spectral range of the reflectivity 
data which makes the Kramers-Kronig transformation 
less reliable. However the use of ellipsometric data ap- 
pears to provide highly accurate results, which are fully 
consistent with the existing EELS data [t2l|]22j |. I(u>) de- 
termined for a number of high-T c superconductors show 
a single broad peak in the energy range below 2 eV. A 
characteristic measure for the width of the peak, which 
can be easily compared with experimental data, is the 
half- width T taken at the low-energy side of the peak at 
xq. For A = 1.0 we obtain T/xq ~ 0.31 (see inset of Fig. 
2). Typical values extracted from experimental loss func- 
tion data are 0.35 for BSCCO, 0.40 for YBa 2 Cu 3 7 , and 
0.45 for Tl2Ba2Ca2Cu30io, i.e. consistent with the the- 
oretical values for T/xq. Furthermore it has been noted 
that the low-energy part of the loss function of all systems 
is quadratic in frequency, i.e. I{u>) oc uj 2 , nearly up to 
the plasma frequency. It has been conjectured that this 
quadratic frequency dependence is universally obeyed in 
all high-T c superconductors |2l],|2^], although at lowest 
frequencies (uj < 0.2 eV) the dependence seems to be lin- 
ear j23|. We note that I(ui < to p ), as given by Eq.(15), 
shows indeed an approximative quadratic dependence, 
which is due to the anomalous damping, Eq.(|l0|), while 
at low frequency the w-dependence becomes linear, be- 
cause of the T-scale in Eq. (|l(i|). 

As we have assumed equidistant planes, we shall focus 
here on the single-layer material LSCO for a quantitative 
comparison between theory and experiment. For moder- 
ate doping the plasmon peak in the loss function is found 
at 0.80 (0.83) eV for x = 0.10 (0.20), respectively [|. 
These values and T/xq ~ 0.3 are in reasonable agreement 
with our numerical results. It is surprising, however, that 
the experimental plasma frequency is not monotonous in 
its doping dependence. For x = 0.34 Uchida et al. || 
report ui p = 0.77 eV. This discrepancy may be related 
to the interplay with interband excitations in the highly 
doped material, which are not contained in the t-J model. 
Yet it is remarkable that in a model where the physics 
scales with the hole-concentration (which is rather small) 
also the correct scale of the plasma oscillation can be ob- 
tained Q. 

In summary we have shown, that the anomalous damp- 
ing of the optical conductivity, as obtained within the t-J 
model, explains the universal line shape of the loss func- 
tion observed in the optimal doping range pi"[ . In this 
regime the width of the plasmon peak is of the order of 



the plasmon energy itself. Model results give besides the 
qualitative also a good quantitative description of EEL 
and the plasmon resonance in cuprates, at least in the 
intermediate doping regime. 
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